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Results from the nuclear recoil calibration of the XENONIOO dark matter detector installed un- 
derground at the Laboratori Nazionali del Gran Sasso (LNGS), Italy are presented. Data from 
• measurements with an external AmBe neutron source are compared with a detailed Monte Carlo 

\l simulation which is used to extract the energy dependent charge-yield Qy and relative scintillation 

^— ' efficiency CeR- A very good level of absolute spectral matching is achieved in both observable signal 

^^ channels - scintillation SI and ionization S2 - along with agreement in the 2-dimensional particle 

\ \ discrimination space. The results confirm the validity of the derived signal acceptance in earlier 

^ reported dark matter searches of the XENONIOO experiment. 

^ PACS numbers: 95.35.-|-d, 14.80. Ly, 29.40.-n, 

5-H Keywords: Dark Matter, WIMPs, Direct Detection, Xenon 

I. INTRODUCTION teracting Massive Particles (WIMPs) off xenon target nu- 

clei, or (in the absence of signal) to set limits on the 
WIMP-nucleon interaction cross-section. XENONIOO is 
The XENONIOO detector [1] aims to detect Galactic a two-phase (liquid/gas) time projection chamber (TPC) 
dark matter through the elastic scattering of Weakly In- with an active volume containing 62 kg of ultra-pure liq- 
uid xenon (LXe), shielded by an active LXe scintillator 

veto containing 99 kg of the same quality liquid. A total 

of 242 1 square inch Hamamatsu R8520 photomultiplier 

♦Electronic address: paul.scovell@physics.ox.ac.uk; present ad- tubes (PMTs) are USed tO read OUt the twO LXe volumes. 

dress: Dept of Physics, Univ. of Oxford, UK The TPC and the active veto are mounted in a low back- 

tElectronic address: marc.weber@mpi-hd.mpg.do ground stainless steel cryostat, enclosed by a passive ra- 

tnot part of the XENONIOO collaboration 



diation shield that effectively attenuates and moderates 
external 7-ray and neutron background [1]. The experi- 
ment is located underground at the Laboratori Nazionali 
del Gran Sasso (LNGS), Italy. XENONIOO simultane- 
ously measures the direct scintillation (SI) and the ion- 
ization signal, via electroluminescence in the gas (S2), 
following recoil energy depositions within its active vol- 
ume. This allows event-by-event particle discrimination, 
and 3-D event vertex reconstruction for definition of a 
fiducial volume and rejection of events from the outer re- 
gions of the TPC. Results from 225 live-days exposure 
have set the most stringent limit on dark matter interac- 
tions to date, excluding spin-independent WIMP-nucleon 
cross-sections above 2 x 10"^^ cm^ for a 55 GeV/c^ 
WIMP mass at 90% confidence level [2]. XENONIOO also 
excludes spin-dependent WIMP-neutron cross-sections 
above 3.5 x 10"''" cm^ for a 45 GeV/c^ WIMP mass [3]. 

This article presents results from the comparison be- 
tween Monte Carlo (MC) simulation of the neutron 
exposure and the equivalent calibration data. Scin- 
tillation and ionization channels are assessed indepen- 
dently as well as simultaneously to compare SI, S2 and 
logio(S2/Sl), the discrimination parameter. This assess- 
ment allows the parameterization of both the charge yield 
Qy, the number of ionization electrons produced by a 
neutron recoil of a given energy, and the relative scin- 
tillation efficiency >Ccffi the energy dependent yield for 
scintillation photons emitted following a nuclear recoil 
interaction. The two quantities are needed to convert 
from observed SI and S2 signals to true recoil energy 
deposition, respectively. The ionization and scintillation 
processes and the correlation between the two are dis- 
cussed in Ref. [4]. 

Finally, given the derived £oS and Qy, studies of ex- 
pected WIMP signals in XENONIOO for two WIMP 
masses at fixed cross-sections are performed. 



A. Energy Calibration 

The energy scale calibration for electronic recoil in- 
teractions in XENONIOO and LXe TPCs in general is 
determined primarily using known 7-ray emission lines 
from standard calibration sources inserted close to the 
LXe volume. Energy deposited by the 7-rays (through 
electronic recoils) leads to characteristic features in the 
physical observable of photoelectrons (PE) in both the 
SI and S2 channels from which the electronic recoil en- 
ergy scale can be inferred. This scale is referred to as the 
electron-equivalent energy scale (keVcc)- 

In the case of nuclear recoils, the scintillation yield is 
quenched with respect to that for electron recoils as typ- 
ically measured at the 122 kcVce 7-ray line of a ^^Co 
source [5]. This can be parameterized by the relative 
scintillation efficiency £off- The observed SI must be cor- 
rected for any spatially dependent effects (such as varia- 
tions in light collection efficiency) as detailed in Ref. [6]. 
This corrected value is given the nomenclature cSl. The 



nuclear recoil energy E (in keVm) is then related to the 
corrected SI by: 

E=S^ V^ — . (1) 

Ly CoS {E) Sni 

Sec is the electric field suppression factor for electronic 
recoils, is measured using the same 122 keVee line and is 
assumed to be energy independent, ^nr is the same quan- 
tity for nuclear recoils and is also assumed to be energy 
independent. In the XENONIOO detector, operating with 
a drift field of 0.53 kV cm^^, ^n^ = 0.95 and S'ee = 0.58 [7]. 
Ly is the light yield at zero field for 122 keVcc electron 
recoils and is determined to be (2.28±0.04) PE/keVco for 
the 225 live-days dark matter search of XENONIOO [2]. 
The S2 signal originates from the ionization electrons 
produced by particles or radiation interacting in the LXe 
volume. The number of free electrons per unit energy is 
the charge yield Qy. The electrons drift towards the liq- 
uid/gas interface in the presence of a homogeneous elec- 
tric field of 0.53 kV cm"^. The ionization signal is de- 
tected as proportional scintillation produced by the ex- 
traction and acceleration of electrons in the gas phase 
above the liquid target volume. A strong electric field 
of ~12 kV cm"^ gives an extraction efficiency of close 
to 100% for an electron reaching the liquid level [8, 9]. 
Thereby, the conversion from extracted number of elec- 
trons to detected S2 photoelectrons is parameterized by 
the secondary amplification factor, Y. As with SI, the 
S2 signal must be corrected for all spatial effects (such 
as electron absorption due to the finite electron lifetime 
in LXe) and is given the nomenclature cS2. Summariz- 
ing both effects, the nuclear recoil energy E (in keVm) 
is related to the ionization signal cS2 (in PE) via the 
following: 



cS2 
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For XENONIOO, Y has been measured by a dedicated 
analysis of the single electron S2 gain. The ampli- 
fied signal follows a gaussian distribution with a well 
defined mean of (19.5 ± 0.1) PE/e" and la width of 
ay = 6.7 PE/e", valid for the time period of the pre- 
sented neutron calibration. A detailed study of the sin- 
gle electron response will be presented in an upcoming 
paper. The conversion between S2 and keVnr may be 
field-dependent (seen in Ref. [7]) but sufficient measure- 
ments have not been made to fully parameterize this de- 
pendence. 

The sensitivity of any dark matter detector depends 
critically on its response to low energy elastic nuclear 
recoil interactions. Typically, broad spectrum neutron 
sources, such as a ^ AmBe (a,n) neutron source, are 
used for calibrating the detector to such nuclear recoils, 
as in XENONIOO. A method for determining the energy- 
dependent response of the detector to nuclear recoils is to 
compare neutron calibration data for a ^ AmBe expo- 
sure with a detailed MC simulation. The ^'*^AmBe cali- 
bration data were taken at the beginning of the reported 



225 live-days dark matter run. This method comple- 
ments direct measurements using dedicated experimental 
setups with lower mass but higher light and charge yields 
per unit recoil energy as presented in e.g. Refs. [10, 11]. 
Historically, the SI and S2 channels have been pre- 
sented independently [12, 13]. However, exploiting both 
SI and S2 channels together in data-MC comparison al- 
lows considerably more robust constraints. In such a way, 
the detector response can be mapped both in cSl, ac- 
cording to the extracted >Ccff, and cS2, with Qy, simulta- 
neously. Such necessary consistency probes the ability of 
the MC simulations to reproduce energy dependent event 
distributions where SI and S2 channels are combined to 
provide discrimination. Additional verification of the de- 
tector response can be achieved through comparison of 
the simulated source neutron emission rate required for 
spectral matching. The neutron emission rate used in 
the simulation must match the independently measured 
source strength for true agreement between data and MC. 



II. MODELING NEUTRON INTERACTIONS 

The XENONIOO instrument including the shielding 
and surrounding environment has been modeled in de- 
tail using the GEANT4 toolkit [14] as previously de- 
scribed in Ref. [15]. The physics inputs for this model 
have been extended to simulate the ^'^^ AmBe neutron cal- 
ibration exposure of XENONIOO, with nuclear recoil an- 
gular cross-sections calculated using the ENDF/B-VI/B- 
VII databases [16] provided in the data files G4NDL 3.13. 
The input ^^^AmBe spectrum adopted is that of an ISO 
8529-1 standard [17]. Results are confirmed to be robust 
against uncertainties in the initial spectral shape since 
the final recoil spectrum in the active volume depends 
only weakly on it. Average variations of less than 5% 
in the simulated recoil energy spectrum are found and 
considered sub-dominant to other sources of uncertainty 
discussed in this publication. 

The neutron emission rate of the ^'^^AmBe source was 
measured at the Physikalisch-Technische Bundesanstalt 
(PTB), the German National Metrology Institute, in Au- 
gust 2012. To measure with high accuracy, the main com- 
ponent of the PTB Bonner sphere neutron spectrometer 
was used. The setup is well suited for (n,7) discrimi- 
nation. The neutron emission rate of the XENONIOO 
^*^AmBe source was determined using the ratio of count 
rates of this source and a reference ^''^AmBe source, 
the rate of which is well known and traceable to na- 
tional standards. The measurement includes a system- 
atic check of the flux isotropy and results in an integral 
source strength of (160±4) n/s. For the purposes of spec- 
tral matching performed in this publication, the source 
strength is fixed by this measured value and its associated 
error and no additional scaling of the simulated neutron 
rate is applied. 

Coincidences between nuclear recoils in the active LXe 
volume with direct 7-rays from the source (such as the 



3.2 MeV or 4.4 MeV 7-rays from de-excitation of -^^C* 
populated by the Be(a,n) reaction, as well as a number of 
low energy 7-rays that do not reach the sensitive volume 
of the detector) are negligible as there is no angular cor- 
relation between the emitted neutron and the coincident 
7-ray. In addition, 7-rays from the ^*^AmBe are attenu- 
ated through the use of a 10 cm thick lead wall mounted 
between the source delivery tube and the LXe volume. 

With the quoted source strength, pile-up effects from 
further source related or delayed emission, such as 
2.2 MeV 7-rays following radiative capture of thermal 
neutrons on the hydrogenous shielding, are also insignif- 
icant (according to the simulation) . Excitation lines due 
to the activation of the xenon by the neutron source give 
signals at high values of logio(cS2/cSl). These are ef- 
fectively removed by cutting on this discrimination pa- 
rameter. In addition, the contribution from scattering to 
metastable states has been studied (as this may lead to 
changes in the spectral state) but has been found to be 
negligible. 

Energy depositions from nuclear recoil interactions in 
the simulation are selected and the process of signal 
generation, including detector threshold and resolution 
effects, is modeled as described in the following sub- 
sections. The final selection of single-scatter recoil events 
includes any multiple scatters that are indistinguishable 
from single scatters due to proximity of vertices in all 
three dimensions, or due to small secondary energy de- 
positions that would be sub-threshold in XENONIOO [1], 
and are present at the level of about 1%. 

Since CcS has been measured with greater accuracy 
and down to lower recoil energies in comparison to 
Qy [10, 11, 13], the following approach to model the sim- 
ulated signal response is adopted: As described below, 
Qy is deduced by performing a y^-^t of the MC gener- 
ated cS2 spectrum to data. The impact of systematic 
uncertainties in correlated parameters, such as signal ef- 
ficiency and £offi are quantified. For the determination 
of Qy, a global fit to all previous direct measurements 
of £off(as described in Ref [18]) is used. Once Qy has 
been obtained, a parameterization of £off is extracted 
using the same y^ minimization technique giving an ab- 
solute matching between the cSl data and corresponding 
MC spectrum. A robust description of both the charge 
and scintillation yield is achieved and, therefore, it is 
shown that the measured cSl and cS2 distributions can 
be matched absolutely to the MC framework using the 
same inputs and methods as developed for recent dark 
matter analyses [6, 18]. 



A. Monte Carlo S2 and SI Signal Generation 

The signal conversion and detector response is com- 
puted based on the energy deposition of nuclear recoil 
events recorded by the GEANT4 model. 

As a first step in generating the eventual cS2 signals, 
the number of electrons n„ liberated at the event vertex 



position is simulated. The mean value fie = Qy {E)E 
is calculated from the definition of the charge yield 
function Qy. Poisson fluctuations of this quantity are 
taken into account on an event-by-event basis. In the 
subsequent step it is considered that the finite electron 
lifetime in LXe [19] determines the number of electrons 
that reach the liquid/gas interface after drifting along a 
uniform electric field. Experimentally, the electron loss 
can be quantified using the depth-dependent charge-yield 
distribution of S2 signals. For the ^^^AmBe calibration 
run, the electron lifetime was measured using the 40 
keVoo signal associated with the de-excitation of ^^^Xe 
nuclei following inelastic interactions with neutrons. 
The exponential suppression factor was found to be 
Tf, = 356 /^s which is consistent with results using a 
"'^^''Cs source (used to determine the electron lifetime 
throughout the dark matter search run). Both methods 
combined have a systematic uncertainty of ~2%. The 
simulation takes into account the signal loss as a function 
of electron drift time td, hence, ne{td) = neexp{-td/Te). 
Finally, the measured Gaussian S2 response of single 
electrons is applied, described by the mean amplification 
Y and width cry, both introduced in the context of 
Eq. (2), to the number ne(td) of simulated electrons 
reaching the gas phase. The resulting signal represents 
the uncorrected S2 and, hence, has all spatial corrections 
removed in order to provide a comparison with the 
raw S2 (S2) recorded in XENONIOO. The procedure 
is repeated in the simulation for every nuclear recoil 
energy deposit recorded along a given neutron track. For 
any two neighbouring scatters at a distance below the 
spatial S2 peak separation resolution of 3 mm along the 
drift direction [I, 18], the S2 signals computed in this 
way are summed because the individual contributions 
are indistinguishable in the real data reconstruction. 
The resulting array of uncorrected S2 signals is ordered 
by size, and for each entry the corrected cS2 value is 
additionally computed by inverting the electron lifetime 
suppression, thus cS2 = S2 exp{td/Te). 

For the generation of SI, simulated recoil energy is 
converted to an observable photoelectron signal using 
Eq. (I). The resolution on this average quantity is 
sampled by assuming a Poisson distribution, which yields 
an integer number. However, since the light collection 
efhciency (LCE) throughout the XENONIOO hquid 
volume is not uniform, simulated numbers of primary 
photons are adjusted using the position dependent, 
experimentally known LCE correction map [I] in order 
to reproduce the true position-dependent light-yield 
and derive the detected, i.e. uncorrected, SI signal 
size. Subsequently, the response of the PMT electron 
amplification to the emission of a single photoelectron 
from the photocathode is taken into account. The mean 
signal resolution of the PMTs, determined using an 
LED calibration to monitor the single-photoelectron 
response and measured to be 0.5 PE for the I PE peak, 
is applied to add Gaussian smearing to the number of 



photoelectrons. Finally, the corrected cSI is inferred by 
applying the experimental LCE map in inverse direction, 
this time inserting the reconstructed event position of 
the largest S2 signal determined before. In this way the 
simulation method correctly resembles the processing of 
actual data. 



B. Deflning Cuts and Efficiencies in Simulated 
Data 



In order to reliably extract Qy and £eff , fiducial, energy 
and data quality cuts are applied to both data and MC 
signals. A fiducial volume of 34 kg (equivalent to that 
used in Ref. [2]) is selected using both radial and drift- 
time cuts. 

It is important to measure the response of XENONIOO 
to single-scatter nuclear recoil events as potential WIMPs 
interacting in the active volume of the detector will do so 
only once. In addition, the SI light from multiple scatter 
interactions will, in fact, result in a single pulse. In the 
multiple-scatter case, the relation between C^s a-nd recoil 
energy in Eq. (I) is inappropriate. Consequently, single 
scatter events are selected in both data and MC using 
identical parameterization as described in Ref. [6]. 

Furthermore, a two-fold PMT coincidence requirement 
for SI signals is defined in the data. Two signals larger 
than 0.3 PE must be seen in separate photomultiplier 
tubes within a time window of ±20 ns. As photon 
tracking is not included in the simulation, the energy- 
dependent efficiency of this cut ~ as determined using 
calibration data [2, 6] - is considered in the MC spec- 
trum. In addition, to avoid threshold effects for small S2 
signals, an S2 threshold (defined using the uncorrected 
electron lifetime S2) of 150 PE is applied to both data 
and MC. 
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FIG. 1: Cut efficiencies used in this analysis. The red 
(dashed) line represents the effect of the S2 threshold effi- 
ciency on cSl. This efficiency is extracted directly from the 
MC simulation, taking the best fit Qy as input. The blue line 
represents the efficiency curve when all data quality cuts have 
been applied. Details of the cut efficiencies used can be found 
in Ref. [6] . In the process of spectral matching, the function 
is allowed to vary within 10%. 



Finally, cuts must be applied to the calibration data to 
remove spurious events that are accepted as single scat- 
ters. No additional noise signals are added to the MC 
simulation, hence, the efficiency of these cuts as derived 
using calibration data is applied to the MC spectrum. 
The definition and energy dependent efficiencies of these 
cuts are discussed in depth in Ref. [6]. 

Fig. 1 shows the efficiency for the S2 threshold cut 
which is extracted directly from the simulation and trans- 
lated to an efficiency as a function of cSl. Also shown 
is the overall efficiency function used in this publication 
which includes all other cuts mentioned above. 



III. METHOD AND RESULTS 
A. Ionization Channel — Determining Qy 

As a ffist step Qy is derived by fitting the simulated 
cS2 spectrum to the one observed in data. In this process, 
CeS remains fixed to the parameterization presented in 
Ref. [18]. 

A x'^-niinimization technique [20] is used to find the 
best matching between data and MC by varying pivot 
points of an Akima spline [21] interpolation of Qy. For 
every intermediate x^ computation, the non-linear de- 
scent algorithm requires the re-evaluation of the detector 
response, applying the updated Qy to generate S2. 

Qy is parameterized by 8 unconstrained and indepen- 
dent spline pivot-points at 0.5, 3, 8, 15, 25, 40, 100 and 
250 keVni-. The lowest pivot point is added to provide 
an unbiased extrapolation to zero recoil energy but has 
effectively no impact on the spectral matching. In data, 
the corrected cS2 spectrum ranges from to 8000 PE, 
divided into 65 bins of equal width. 

The impact of various simulation parameters on the 
best-fit Qy was studied to estimate the systematic un- 
certainty of the final result. The largest systematic error 
is connected to the choice of CcS as variations in this 
quantity lead to changes in the simulated cSl spectrum 
and, consequently, in the number of events passing the 
selection requirements. With a lower (higher) value of 
£eff the cSl energy spectrum of accepted events will be 
shifted upwards (downwards). Accordingly, Qy will de- 
crease (increase) in order to compensate this effect and 
re-establish the matching in cS2. This interdependency 
is present mainly near the detection threshold, where the 
acceptance as function of cSl falls steeply (Fig. 1), and 
becomes negligible at higher recoil energies. The Ccs pa- 
rameterization is allowed to vary within the ±1ct uncer- 
tainty bounds as defined in Ref. [18]. Similarly, the cSl 
efficiency function was allowed to vary by ±10% around 
its reported mean. The systematic error connected to 
the choice of pivot positions and initial values has been 
found to be negligible in the energy region above 3 keVnr 
(the lowest energy at which £off has been directly mea- 
sured [10]). Finally, the statistical uncertainty of about 
1% on average is also included. This is obtained after 



repeating the simulations about 50 times at fixed config 
urations but varying random seeds. 
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FIG. 2: Comparison between the MC and data cS2 spec- 
tra. The black data-points indicate the data and the blue 
spectrum is obtained as the result of the optimization of Qy. 
Good agreement between spectral shape and absolute rate 
across the whole signal range is achieved. For comparison, 
the gray dashed line indicates a generated cS2 spectrum, as- 
suming the same Qy as shown by the dashed line in Fig. 3 
and described in Ref. [22]. 
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FIG. 3: Result on Qy obtained from fitting the MC generated 
cS2 spectrum to data. Pivot points of the spline interpola- 
tion are shown in light blue. The shaded area indicates the 
systematic uncertainty from varying input parameters of the 
simulation (find discussion in text). The interpolation be- 
tween the pivot points at 0.5 and 3 keVnr does not yield a 
reliable result for Qy and is shown using a dashed-blue line. 
The purple data points show the result of the first measure- 
ment of Qy in LXe at 0.2 kV cm"^ [7]. Red data points 
show the result from direct measurements at a drift field of 
1.0 kV cm"^ [11]. The green hatched area is the combined 
result from the ZEPLIN-III experiment, extracted in a simi- 
lar fashion to this work although at a much higher field [13]. 
The black dashed line represents a predicted Qy based on a 
specific phenomenological model as described in Ref. [22]. 

The resulting pivot points and systematic errors to- 
gether with the spline interpolation yield a best-fit Qy 
function. Fig. 2 shows the spectral matching correspond- 
ing to the central fit value of Qy (shown in Fig. 3) along 



with data. The spectral matching is good from low to 
high values of cS2. In Fig. 3 Qy is determined down to 
~3 keVnr with similar precision as achieved in the direct 
measurement [11] but, in some energy intervals, devia- 
tions of the central value from the direct measurement of 
Manzur et al. are observed. Best matching is achieved 
with a simulated neutron emission rate of 159 n/s which 
is in agreement with the independently measured neutron 
emission rate of (160 ± 4) n/s as described in Sec. II. 



B. Scintillation Channel — Determining jCcH 

As with the extraction of the best-fit Qy in the previ- 
ous sub-section, the £eS fit is parameterized by the same 
8 independent spline points, which are allowed to vary to 
give the best agreement between data and simulated cSl 
spectrum. The fit range is constrained to 2-200 PE as 
good agreement below 2 PE has not been achieved for a 
wide variety of £eff parameterizations. This mismatch is 
predominantly due to uncertainties in the calculated effi- 
ciency curve as given in Ref. [2] but it may, in part, also 
be due to uncertainties in the neutron physics provided 
by GEANT4 at the lowest recoil energies. A further cause 
of uncertainty is that the calibration of the detector re- 
sponse becomes difficult as signals approach the single- 
or two-photoelectron level where the PMT response due 
to single photoelectron size noise or electronics noise be- 
comes more difficult to characterize. 

For the simultaneous generation of S2 as a part of the 
complete signal simulation, the central curve of the previ- 
ous fit result for Qy is used. Consequently, any result on 
the best-fit Ccs has a systematic uncertainty of the size 
of the ±1(7 bounds of the CcS representation in Ref. [18], 
which was found to be the main contributor for the esti- 
mated error on Qy as presented in Fig. 3. 

Fig. 4 presents the result of matching for the £eff (gray) 
as described in Ref. [18] and applied at the beginning 
of the fitting process while the best-fit £eff line is rep- 
resented by the blue spectrum. The optimized Ccs is 
shown in Fig. 5 in comparison with the literature values. 

It is shown that this CeS is in good agreement with 
that measured by Plante et al. [10] and in good agree- 
ment overall with other measurements below 15 keVnr- 
The deviation between the extracted CcS and the mean 
measurement results from the improvement of the spec- 
tral matching in the range of 20 - 60 PE in Fig. 4. 

The mean extracted £off provides an important con- 
sistency check to recent direct measurements but the 
method does not improve on the accuracy on this quan- 
tity until more precise direct measurements of Qy as a 
function of energy and drift field become available. 



C. Two-Dimensional Distributions 

When satisfactory data/simulation agreement is 
achieved for the ionization and scintillation channels in- 
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FIG. 4: Fit of the simulated cSl spectrum to data (black 
points). The MC spectrum (blue) is obtained using the £cff 
after the optimization process. Below 2 PE (indicated by the 
vertical red line), a discrepancy between data and MC is ob- 
served (see text for discussion). The gray dashed Hne shows 
the spectral shape using the L^g detailed in Ref. [18] for com- 
parison. Reasonable agreement between data and MC above 
2 PE is already achieved with this jCcH parameterization. 




FIG. 5: yCeff (blue line) obtained after the optimization of 
the absolute cSl matching. As with Fig. 3, pivot points of the 
spline interpolation are shown in light blue. As with the ex- 
traction of Qy, the parameterization of Leff is unreliable below 
3 keVnrand is shown using a dashed blue line. For compari- 
son, literature values of jC^b including Aprile et al. (•) [23], 
Manzur et al. (A) [11], Plante et al. (■) [10], Horn et al. 
combined result (green-shaded) [13] are shown along with this 
work (blue). Also shown is a global fit to all LeB data, used 
in Ref. [18] (black line and gray-shaded uncertainty). 



dividually, the combined two-dimensional distributions 
provide a robust test of the consistency of the derived 
nuclear recoil energy scales. This is because even though 
individually the cSl spectrum is only weakly sensitive to 
changes in Qy and the cS2 scale only weakly sensitive to 
changes in £cff j a two-dimensional comparison is sensitive 
to both. 

The logio(cS2/cSl) projections of the two-dimensional 
distributions for both data and MC are sliced into 2 PE 
bins and compared. Fig. 6 shows several of these distri- 
butions. Upon examination, it is clear that the match- 




(such as between PMTs in the lower array) and cannot 
be predicted. 
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FIG. 7: Two-dimensional distributions of the means of the 
single scatter elastic nuclear recoil population from AmBe 
calibration data (black) and MC simulations (blue), where 
log(cS2/cSl) is the discrimination parameter and the energy 
scale is defined by the primary scintillation channel (top). 
The error bars (black) and filled area (blue) represent the 
±1(7 bands for data and MC respectively. The bottom panel 
shows the residual differences between data and MC which 
are all within 2%. 



FIG. 6: Projection of logio(cS2/cSl) in data (black) and 
MC (blue) for slices in SI. Good agreement is found for 
all slices. Mild deviations are observed at lower values of 
logio(cS2/cSl). A possible reason for these deviations is dis- 
cussed in the text. 



ing between data and MC is reasonable but there are 
some variations, particularly at the lowest values of 
logio(cS2/cSl). This variation is most likely due to the 
presence of anomalous events containing two recoil ver- 
tices (cSl and cSl'), one of which (cSl') occurs in a re- 
gion where the electric field configuration will not allow 
ionization electrons to be drifted to the liquid/gas in- 
terface. The recoil which occurs in such a region will 
have no associated ionization signal and, as such, is in- 
distinguishable from a single scatter event (for electronic 
recoil backgrounds, this effect has been reported initially 
in Ref. [24] and subsequently in Ref. [25]). These sig- 
nals will have a relatively larger cSl signal as compared 
to that expected given the size of the cS2 signal (actu- 
ally associated to only one of the two recoils) causing 
the apparent ratio of cS2/cSl to fall. The simulation 
includes events where the second recoil occurs between 
the cathode and the lower PMT array. The cSl' signal 
is calculated using a LCE map which has been extrapo- 
lated below the cathode. Any remaining variations can 
be attributed to recoil events where the cSl' recoil oc- 
curs in a region where the LCE is not precisely known 



Fig. 7 shows the parameterization of the elastic nu- 
clear recoil band from the ^''^AmBe calibration in data 
and MC in the phase-space of logio(cS2/cSl) against cSl. 
Points represent the medians of the distribution intervals 
in energy as defined by the scintillation channel. Vertical 
bars represent the ±la quantiles of each slice, encom- 
passing the statistical fluctuations in signal generation 
and detector resolution. Good agreement is achieved 
between calibration data and MC for both means and 
widths of each distribution. Although the logio(cS2/cSl) 
distribution is systematically higher for the MC, it is 
still within 2% of the central data values. In this 2- 
dimensional distribution, the effect seen at low values 
of logio(cS2/cSl) in Fig. 6 will cause events to move to 
logio(cS2/(cSl-|-cSl')) whilst simultaneously moving to 
a higher apparent value of cSl (namely cSH-cSl'). The 
most pronounced changes will occur when cSl' is of a 
similar or greater size than cSl. The level of agreement 
between data and MC in this comparison clearly shows 
that the effect of multiple scatter events that appear con- 
sistent with a single scatter event is minimal (as the max- 
imum observed deviation is < 2%). 

In addition to comparing logio(cS2/cSl) vs cSl for 
data and MC, a comparison in cS2 vs cSl is shown in 
Fig. 8. A good level of agreement is demonstrated over 
the full range in cSl and cS2. Disagreement between 
data and MC increases at higher energies as the recoil 
energy spectrum from ^ AmBe falls exponentially with 
recoil energy. 
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FIG. 8: Comparison of cS2 vs cSl for data (black) and MC 
(blue). The solid lines represent the median of the cS2 vs 
cSl distribution while the vertical lines (black) and filled area 
(blue) represent the ±la bands for data and MC respectively. 
As with the parameterization reported in [18], the extracted 
£cs infers a energy scale that reaches 1 keVnr at PE. 



The level of matching achieved in the 2-diniensional 
comparison confirms that it is possible to reliably predict 
nuclear recoil event signatures in the both signal channels 
simultaneously. 



IV. SIMULATED WIMP DISTRIBUTION 

Simulating the response of XENONIOO in both ioniza- 
tion and scintillation channels gives the ability to predict 
the distribution of WIMP recoil events in the discrimi- 
nation parameter space. An isothermal WIMP halo with 
a local density of p^ = 0.3 GeV cm"'^, a local circular 
velocity of i/q = 220 km s"^, and a Galactic escape veloc- 
ity of 544 km s" [26] is assumed. For the calculation of 
expected rates, CcS from Ref. [18] and Qy as determined 
in this work are used. 

Sample recoil spectra for both an 8 GeV c"^ WIMP 
mass with a spin-independent WIMP-nucleon cross- 
section of 3 X 10-41 cm^ and for a 25 GeV c'^ WIMP mass 
with a spin-independent WIMP-nucleon cross-section of 
1.6 X 10"'''^ cm^ are shown in Fig. 9. WIMP masses are 
chosen such that the signal response to both light and 
more massive particles is investigated. The choice of the 
interaction cross-sections is motivated by the intersec- 
tion of these masses with the lower boundary (central 
region) of recent dark matter signal claims from the Co- 
GeNT [27] and CRESST-II [28] experiments. The distri- 
butions are generated for a dark matter exposure equiva- 
lent to that of the 225 live-days XENONIOO dark matter 
search. From the given recoil spectra, simulated cSl and 
cS2 signals are generated whilst maintaining the simu- 
lated cuts described in Sec. II B. The efficiency for all 
data quality cuts is also applied. 

Fig. 10 shows simulated distributions in discrimination 
parameter space for the two WIMP masses superimposed 
on the XENONIOO 225 live-days dark matter search data 
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FIG. 9: Recoil spectra for (solid blue) an 8 GeV c ^ 
WIMP and (dashed blue) a 25 GeV c"^ WIMP with spin- 
independent WIMP-nucleon cross-sections of 3 x 10"*^ cm^ 
and of 1.6 x 10" cm , respectively, as they would be ob- 
served in XENONIOO taking into account energy resolution 
and detection efficiencies. Boundaries of 3 and 30 PE in cSl 
are marked by the dashed red lines. In both cases, standard 
WIMP parameters are assumed (see Sec. IV) and an exposure 
equivalent to that of the 225 live-days XENONIOO WIMP 
search run with a 34 kg fiducial mass is used. 



for comparison [2]. The distributions are presented in a 
parameter space which is flattened by subtracting the 
mean of the electronic recoil distribution. As expected, 
the recoil spectrum of a 25 GeV c"^ WIMP gives rise to 
a larger number of events extending to higher recoil ener- 
gies than the 8 GeV c"^ WIMP and its bulk distribution 
resembles very well the shape of broad-band ^"^^AmBe 
neutron calibration. For the lower mass WIMP, the ma- 
jority of events are expected below the analysis thresh- 
old in a region where XENONIOO is still sensitive, but 
a significant number of recoil events fall above the im- 
posed analysis threshold of 3 PE. However, the center of 
the signal distribution is clearly shifted towards a lower 
login (cS2/cSl) with respect to the average neutron band 
position. This can be explained by the steeply falling 
recoil energy spectrum of light mass WIMPs combined 
with asymmetric upward fluctuations in the Poisson dom- 
inated regime of small SI. For both cases the expected 
number of WIMPs is calculated for an exposure equiva- 
lent to the 225 live-days XENONIOO WIMP search run 
in a region given by a SI range of 3-30 PE SI and below 
the 99.75% electron recoil line as defined in Ref. [2]. The 
results are 270^1^° (sys.) and 1590 ± 40 (stat.) events 
for the 8 GeV c'^ and for the 25 GeV c'^ WIMP, re- 
spectively. For the 25 GeV c"^ case, the statistical error 
exceeds all other systematic uncertainties whereas in the 
8 GeV c"^ case, the expected rate is dominated by the 
systematic uncertainties in e.g. C^s arid Qy for lower en- 
ergy recoils. Rates could similarly be calculated for the 
CcS extracted in this publication. The shape of this Leff 
leads to predicted rates consistent with those calculated 
for the direct Leff measurements within errors. 



The excess of predicted WIMP recoil rates above 
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FIG. 10: Two-dimensional distributions of expected cSl and 
cS2 signals for (top) an 8 GeV c"^ WIMP and for (bottom) 
a 25 GeV c"^ WIMP with spin-independent WIMP-nucleon 
cross-sections of 3 x 10" ^ cm and of 1.6 x 10" cm , re- 
spectively. In both cases, the same assumptions applied to 
create the recoil spectra in Fig. 9 are used. In both figures, 
the vertical red lines represent boundaries of 3-30 PE. The 
horizontal (long-dash) red curve represents the mean (fi) -3a 
for the elastic nuclear recoil distribution and the horizontal 
(short-dash) red curve represents the 99.75% electron recoil 
rejection line as discussed in Ref. [2]. 



only 2 event candidates observed in the 225 live-days 
XENONIOO dark matter search [2] is consistent with the 
reported exclusion limit, supporting the tension between 



these results and signal claims by other experiments [27- 
29]. 



V. CONCLUSIONS 

The neutron calibration of the XENONIOO dark mat- 
ter detector with a ^^^AmBe source has been modeled 
with a MC simulation that includes the signal genera- 
tion in both the SI and S2 channels. Agreement in the 
ionization channel is achieved through the adoption of a 
Qy (derived using ^"^^AmBe data and a fixed £off) that is 
largely consistent with previous direct and indirect mea- 
surements and phenomenological estimations but shows 
no indication of a low-energy increase as reported by the 
direct measurement of Ref. [11]. Additionally, an op- 
timized £cff is determined using a similar method and 
is used to match data and MC signal distributions in 
the scintillation channel. The ionization and scintilla- 
tion channels are combined in two-dimensional spaces, 
achieving agreement between MC and data, constraining 
the uncertainty in the nuclear recoil energy scales, and 
reproducing both means and widths of energy distribu- 
tions. It provides a strong validation of the understand- 
ing of the discrimination parameter space in which previ- 
ous XENONIOO dark matter searches were analysed and 
reported. A simulated neutron emission rate of 159 n/s 
is required to achieve spectral matching. This is in agree- 
ment with the measured emission rate of (160 ± 4) n/s 
and confirms the robustness of the SI signal acceptance 
used in the XENONIOO WIMP searches [2, 3, 6, 18]. 
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